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ABSTRACT 


A  theory  of  ferromagnetism  for  general  spin,  approximately 

valid  through  the  entire  temperature  range,  is  given.  At  low 

temperatures  the  magnetization  agrees  with  the  Dyson  results, 

3  4 

having  no  term  in  1r  and  having  a  term  in  T  equal  to  that 

fourd  by  Dyson  in  first  Bom  approximation;  terms  arising 

from  the  approximations  of  the  theory  first  appear  in  order 

^3(2S  +  ly^  that  a  spurious  T^  term  does  appear  for  S  =  1/2, 

but  for  no  other  spin.  Curie  temperatures  are  within  a  few 

percent  of  the  Brown  and  Luttlnger  estimates  for  spins  greater 

than  unity,  and  agree  within  one  percent  of  the  Domb  and 

Sykes  estimate  of  the  rge  spin  limit.  The  susceptibility  at 

high  temperatures  agrees  with  the  Opechowski  expa:ision  to 
1/  2  ~ 

terms  in  'T  .  The  quasi-particle  energies  are  renoimiallzed 
by  the  energy  at  low  temperature  and  by  the  magnetisation 
at  higher  temperature.  The  Green  function  is  decoupled  by  a 
physical  criterion  involving  self-consistency  of  the  decoupling 
at  all  temperatures.  The  Green  function  method  is  extended 
to  higher  spin  by  a  technique  of  parametrizing  the  Green 
function  and  expllcltely  finding  the  functional  dependence 
on  this  parameter  by  solution  of  an  ataxlliary  differential 
equation. 


Introduction 


The  Heisenberg  model  of  a  ferromagnet  has  been  theoretically 

analyzed  by  Dyson by  series  expansion  in  powers  of  T,  valid 

(o\ 

at  low  temperatures,  and  by  Opechowski'  '  by  series  expansion 
1 

in  valid  at  very  high  temperatures.  We  here  develop  an 
approximate  theory  which  covers  the  entire  temperature  range, 
including  the  particularly  interesting  Intermediate  region 
in  the  neighborhood  of  the  Curie  temperature,  and  which  agrees 
satisfactorily  with  the  rigorous  restilts  at  both  very  low 
and  very  high  temperatures. 

At  low  temperatures  the  magnetization  has  terms  of  order 

5/2  5/2  74 

T  ,  T  ,  T  which  agree  with  the  Dyson  results,  it  properly 

3  li 

has  no  term  in  T  ,  and  the  term  in  T  is  equal  to  that  found 

by  Dyson  in  first  Born  approximation,  T3rms  arising  froim  the 

.  3(2S  +  l)/2. 

approximation  in  the  theory  first  arise  in  order  T  , 

so  that  the  case  of  S  =  1/2  is  an  exception  to  the  above 
statements,  having  a  spurious  T  term.  Curie  temperatures 
are  quite  close  for  spin  2)  to  the  values  estim.ated  by 

Brown  and  Luttlnger by  extrapolation  of  the  high  temperatures 
series  expansion,  except  for  very  small  values  of  spin  (again 
the  case  of  spin  1/2  is  particularly  unsatisfactory).  The 
Curie  temperatures  agree  within  1^  with  the  estimate  of  Domb 
and  Sykes^^4or  the  high  spin  limit.  The  susceptibility  at 

(2) 

high  temperatures  for  all  spin  values  agrees  with  the  Opechowski''  ' 


1/  2 

expansion  to  terms  In  ^  T  . 

The  quasi-particle  energies  are  equivalent  to  simple 
spin  wave  energies  "renormalized"  by  a  factor  which  is 
proportional  to  the  thermodynamic  energy  at  low  temperatures, 
but  which  becomes  proportional  to  the  magnetization  at  higher 
temperatures. 

An  heuristic  interpretation  of  the  renormalization  of 
quasi-particle  energies  has  been  given  by  Keffer  and  Loudon ' . 
They  point  out  that,  at  low  temperatures  where  only  long  wave 
length  spin  waves  are  excited,  the  local  magnetization  direction 
varies  slowly  through  the  crystal.  Excitation  of  an  additional 
spin  wave  is  analogous  to  excitation  of  a  ripple  relative  to 
this  slowly  varying  local  magnetization.  The  effective 
exchange  Integral  determining  the  energy  of  this  ripple  is 
influenced  by  the  angle  between  neighboring  spins  in  the 
slowly  varying  background  medium.  This  angle  also  determines 
the  thermodynamic  energy;  hence  the  renormalization  of  the 
spin  wave  energy  by  the  thermodynamic  energy.  However  at 
higher  temperatures  the  thermally  excited  excitations  have 
wave  lengths  comparable  to  the  Inter-spln  distance,  and 
the  correlation  distance  in  the  background  medium  is  as  short, 
or  shorter,  than  the  wave  length  of  the  particular  excitation 
being  considered.  In  this  region  our  results  indicate  that 
the  effective  exchange  integral  is  Influenced  by  the  angle 


relative  to  the  average  magnetization;  that  is,  the  spin  wave 
energies  are  renormalized  by  the  average  magnetization. 

The  Heisenberg  ferromagnet  with  spin  1/2  was  analyzed  by 

Tyablikov^^^  using  the  technique  of  double-time  temperature- 

dependent  Green  functions  Extension  of  the  theory  to 

higher  spin  has  been  achieved  recently  by  Tahir- Kheli  and 

(q)  "■ 

ter  Haar^^'.  The  present  theory  differs  from  those  applica¬ 
tions  of  the  Green  function  method  in  two  respects.  First, 
the  decoupling  of  the  higher  order  Green  functions  is  guided 
by  a  piavisible  physical  criterion.  Second,  the  method  of  employing 
Green  functions  for  general  spin  is  simplified;  the  Green 
function  is  parametrized  and  the  functional  dependence  on 
this  parameter  is  found  explicitely  by  solution  of  an 
auxiliary  differential  equation.  This  single  differential 
equation  replaces  the  2S  simultaneous  Green  function 
equations  in  the  method  of  Tahlr-Khell  and  ter  Haar^^^. 


2.  The  Green  Function 


‘i 


The  problem  to  be  considered  is  the  statistical 
mechanics  of  the  system  described  by  the  Hamiltonian 

W  ®  H  H  )  S.  •  (1) 

^  ^  3.f  ^ 

where  M-S  is  the  magnetic  moment  per  ion;  H  is  the  applied 
magnetic  field  (in  the  positive  z-direction) ;  ^  is  the 
spin  operator  (in  units  of  "fe)  for  the  ion  at  site  g,  and 
j(g-f)  is  the  exchange  j.ntegral  between  ions  at  sites  g  and  f. 
The  exchange  integral  is  assimed  to  be  a  function  only  of  the 
distance  between  ions;  it  is  not  restricted  to  nearest 
neighbors  or  to  only  positive  values,  but  it  is  assumed  that 
the  ground  state  of  the  system  is  a  ferromagnetically  aligned 
state. 

The  temperature-dependent  retarded  Green  function 
involving  the  two  operators  A  and  B,  «A(t);  B»,  is  defined 
by  (8) 


(2) 


Where  A(t)  is  the  Heisenberg  operator  at  time  t;  ^(t)  is 


unity  for  positive  t  and  zero  for  negative  t;  square 

brackets  denote  a  commutator;  and  single  angular  brackets 

denote  an  average  with  respect  to  the  canonical  density 

matrix  of  the  system  at  temperature  T.  The  Fourier  transform 

of  the  Green  fvinctlon  Is  a  function  of  lo  (or  of  E 

and  Is  denoted  by  It  satisfies  the  equation  of  motion 

(3) 


If  this  equation  can  be  solved  for  one  then 

extracts  knowledge  of  the  correlation  fimctlon  '^BA(t)>  from 
the  relation 


^  i!  f  ' P^.zi 


-i  u»t 

(2.  doj 


(M 


Equations  (5)  and  (4)  are  the  only  equations  required  for 
the  application  of  the  Green  function  method. 


For  reasons  which  will  become  evident  subsequently  vre 


consider  the  Green  function 

S  ^  (5) 

where  a  is  a  parameter.  The  Fourier  transform  of  this  Green 
function,  G^  (g,!)^  satisfies  the  equation  of  motion  (eqn.  3) 


(6) 


where 


>^ws<[s",e  5-]> 


(7) 


The  commutator  of  with  the  Hamiltonian,  required 

o 


in  the  last  term  of  equation  (6),  is  easily  computed,  giving 


f  y  ;e s;»^ 


(8) 


The  remaining  problem  is  to  express  the  higher-order  Green 
function  on  the  right  in  terms  of  two-particle  Green  functions, 
so  that  equation  (8)  can  be  explicitely  solved  for  G^  (g,l). 


3.  A  Decoupling  Approximation 

The  (Essential  approximation  in  the  methods  of  Tyabllkov  and  of 
Tahir-Kheli  and  ter  Haar  consists  of  Ignoring  the  fluctuations  of 
S^,  replacing  this  operator  by  Its  average  value: 

s 

(9) 

This  approximation  results  In  the  magnetization  renormalization  of 
quasi-particle  energies,  in  disagreement  with  the  low- temperature 
theory. 

The  decoupling  approximation  to  be  used  here  Is  most  clearly 
described  for  the  special  case  of  spin  In  that  case  we  can 
write  Sg  In  either  of  the  following  forms 


I 


5  5^ 


(S=>/2.) 


s'  -  5  ^ 


(10) 

(ir. 


(5=  Vi')  (12) 


or,  multiplying  the  first  of  these  equations  by  an  arbltrai?y 
parameter  c<  and  the  second  by  (1  -U  )  and  adding, 

C  I  C  ^  C  "  ^  r'  ~  C  ^ 

^  I  \  ^ 

The  Green  function  «5;  s;  5/  '  is  reasonably  de- 

*  o  T  ' 

coupled  in  the  symmetric  form* 

<s,'s;5;;b»  Ks;s;x<ss^•8»  («> 

and  similarly  for  the  Green  function  «s;s:s*\e.». 

0  ^  i 

Thus  the  identity  (12)  leads  to 


*The  remaining  "contraction” 

+  + 

because  the  operator  S  S„  1 

g  f 


vanishes 

s  not  diagonal  In  the  total  z-com- 


ponent  of  spin. 


r 


s;;b»  --  -  «<s;s;>«s;;  b»  w 

Ifo(is  chosen  as  unity  the  result  corresponds  to  decoupling  on  the 

basis  of  identity  (10);(5(.=  o  corresponds  to  decoupling  on  the  basis 

of  identity  (ll);c<=  -1  corresponds  to  decoupling  on  the  basis 

of  the  identity  =  -  S  +  S”*”  S".  Thus  we  are  faced  with  the 

g  g  g 

possibility  of  obtaining  a  correction  to  the  Tyablikov  decoupling 
with  either  a  positive  or  negative  sign,  or  no  correction  at  all, 
or  any  intermediate  value,  depending  on  the  choice  of  o(.  Clearly 
a  physical  criterion  is  required  at  this  point. 

The  operator  S”  in  equation  (10)  represents  the  deviation 

Z  —  “h 

of  S  from  +S.  It  is  this  operator  S  S  which  is  treated  approxi¬ 
mately  when  decoupling  on  the  basis  of  equation  (10).  It  therefore 
seems  reasonable  to  use  equation  (10)  when  the  deviation  from 
=  S  is  small;  that  is, when  ~  S. 

Similarly,  the  operator  ^  (S’^S"  -  S~S^)  in  equation  (11) 
represents  the  deviation  of  S^  from  zero,  and  it  therefore  is 
reasonable  to  use  equation  (11)  as  the  basis  of  decoupling  when 

<  S'>  ^  0  . 


Both  of  the  above  observations  are  contained  in  the  choice 

(15) 


V  = 


(S-'/i) 


for  then  equation  (12)  becomes 

S3  -<S  >  +L  IS  18  3  3  j 


(16) 


The  operator  in  brackets,  which  is  to  be  decoupled 


represents  the 


deviation  of  from  ^S^j>and  should  be  self-consistently  small  in 
all  temperature  regions. 

Inserting  the  above  value  of  (/into  the  decoupling  equation  (19) 
gives 


This  is  the  basic  decoupling  approximation  for  spin  -jj  we  now 
generalize  it  for  higher  spin. 

The  analogue  of  equation  (10),  for  general  spin,  is 


(17) 


=  5(s+i')  ('S^)  -  (18) 

whereas  equation  (11 )  remains  true.  Decoupling  as  before,  and 

z  2 

neglecting*  the  fluctuations  of  (S  )  ,  we  find  in  this  case  as  well 


<(S,* S; ; B»  ^ 8»  (19) 

where «Kls  the  fractional  contribution  of  the  identity  (l8),  and 
(l-o<)  is  the  contribution  of  the  identity  (11 ),  to  this  result. 
Unfortunately  is  no  longer  the  only  operator  treated  ap¬ 

proximately  in  decoupling  equation  (l8),  and  the  Interpretation 
of  the  decoupled  operator  as  being  the  deviation  from  =  +  S 
is  no  longer  true.  Hence  the  choice  of^ls  no  longer  quite  so 
evident.  However, 


*The  Green  function  also  can  be  symmetrically 

decoupled,  but  the  results  are  of  the  same  form  as  those  obtained 


by  the  simpler  procedure  above. 


/6 


we  determine  it  by  the  following  requirements: 

a)  For  S  =  -^  ,  should  reduce  to  our  previous  result^or  to 

-  <S^>/S. 

b)  For  =  o,y  should  vanish.  This  requirements  follows 

from  the  fact  that  Identity  (11)  retains  its  interpretation 
for  arbitrary  S. 

c)  For  ^  S  we  expect  that  should  have  the  form 

z 

S  —  S  -  n,  where  n  is  a  deviation  which  is  of  order  unity 
rather  than  of  order  3. 

Requirement  c)  Implies  that  should  be  of  order  unity, 

rather  than  of  order  3,  at  low  temperatures.  Now  3'*’^  is  the 

‘  s  s 


spin  deviation  in  lowest  order  ;slmllarly  ■<(3"  3^^  will  be  of  order 


23  if  f  and  g  are  closely  coupled, 

1 


Hence  we  take 


o<  = 


(20) 


and  it  is  clear  that  this  satisfies  all  of  the  physical  requirements 
above.  Equations  (19)  and  (20)  characterize  the  basic  decoupling 
approximation  of  the  theory. 


*  The  appearance  of  23  here  corresponds  to  the  identification  of 
3"^  and  3  respectively  as'|23  a"^  and  ^23"  a  in  the  leading'  terms 
both  of  the  Dyson  and  the  Holstein-Prlmakoff  transformations.  Here 
a"^  and  a  are  destruction  and  creation  operators  of  elementary 


boson-type  excitations. 


n 


4 .  Solution  of  the  Green  Function  Equation 

Inserting  the  decoupling  approximation  Into  the  equation  of 
motion  (8)  gives 

These  equations  are  a  set  of  coupled  equations  for  various 
pairs  of  sites  (g,  1),  (f,  l).  Translational  invariance  dictates 
consideration  of  the  spatial  Fourier  transforms 


(22) 


r- 

J(I^)  -1  e  J(y^) 

^  c"  c'^\ 

<e  \ 


(^3) 


(24) 


Here  g*k  denotes  the  vector  product  R *  where  R  is  the  position 

S  S 

vector  of  the  g-site.  Equation  (21)  then  Implies 

(S’)  + 


or 


h<J^  £  [  JriT'i  -0‘(  ic'~x)] 

<' 

hQi  (A) 

'  Tr  (  E-E(i2))' 


(2q) 


(2b) 


where 


K' 

3.  '  ^ 

The  correlation  function  which  can  be  obtained  from  ( k)  by 
equation  (4)  is  uhe  quantity'^(  k,  a ),  defined  in  equation  (24).  In¬ 
serting  equation  (2b)  in  (4),  and  taking  t  =  o,  we  thereby  obtain 


(2ti) 


Equation  (2b),  together  with  eq.  (27)  defining  E(k),  equation 
(24)  def Inlng ^( k> a ) ,  and  equation  (7)  def lnlngt^( a ) ,  is  the  basic 
equation  of  the  theory.  It  must,  however,  be  augmented  by  an  ex¬ 
plicit  relationship  between^ and  *©•,  and  it  is  this  step  which 
complicates  the  problem  for  spln>i. 

The  problem  which  distinguishes  the  simple  case  of  spin  from 
the  more  complicated  case  of  higher  spin  becomes  evident  if  we 
particularize  our  solution  temporarily  to  S  =  ■^.  We  also  take 
a  =  o,  and  we  note  from  eqs.(24)  and  (10)  that 

=<C5''>  =  (S-'Vd 

k  o  3  ^  )  (29) 

Thus^(o)  is  simply  related  to  the  magnetization.  N^(o)  is  the 

— * 

total  number  of  spin  reversals  in  the  crystal,  and'^(k,  o)  is  the 
occupation  number  of  the  basic  excitation  of  wave  vector  K 
Similarly,  from  equation  (7 ) 


(50) 


<£ —  Hence  equations  (28)  and  (50)  determine  ^  (k,  o )  as  a  function 
of  ^  ,  and  equation  (29)  provides  a  requirement  of  self-consistency 
which  determines  . 

For  higher  spin  values  equation  (29)  is  replaced  by 

^(0)1  =<S'S'->  =${S41')-<(S*'|'>  -<s*^ 

This  Introduces  the  new  quantity  <(sb">  ,  which  is  not  determined 

by  equation  (29).  Tahlr-Khell  and  ter  Haar^^ ^therefore  introduce 

the  Green  function  to  evaluate  this  quantity,  but  its 

solution  introduces  <(^(  )^  ^ .  By  introducing  23  such  Green  functions, 

z  2  S”f*l 

and  by  recalling  that  (S  )  is  related  by  an  identity  to  lower 

powers  of  S^,  they  are  thereby  able  to  obtain  a  solution.  That 
method  gets  laborious  very  quickly,  but  they  have  given  explicit 
solutions  to  3  =  3- 

By  exploiting  the  functional  dependence  of  \^a)  and^on  the 
parameter  a,  which  was  inserted  in  equation  (p)  for  just  this 
purpose,  both and V can  be  explicitly  related  to  for  arbitrary 

spin,  in  close  analogy  with  the  case  of  spin 


t).  Relationship  of  #to 

Calculation  of  the  commutator  def Inlng  t6'(a)  (eq.  (7))  Is 
facilitated  by  the  Identity 

[S'*,  (S')"]  -  {(S'-If -(s‘f  i  S»-  (5^) 

which  Is  easily  corroborated  for  n  =  1,  2  and  extended  to  higher  n 


by  mathematical  Induction.  It  follows  that 

^  \  s' 


[S*  -  (e'-Oe^^ 


and  thence 


,  ^<e‘^s'>  +(e"-i)<e*'Vs-> 

Expressing  S’^S“  In  terms  of  by  the  Identity  (l8). 

Finally,  It  Is  convenient  to  Introduce  the  quantity 


AM  =  <e‘^> 


and  the  notation 


D  =  -^ 


da 

Then  the  two  quantities  of  Interest  can  be  written 

‘  S(ui)fe"-i)jl  V  (e'Vi)Dn-(e“'-i) 

and 


(53) 


(54) 


(56) 


(36) 


(5Y) 


(58) 


(59) 


The  relationship  (2ti)  between  V(  a- )  andi6t(a)  can  be  recast 
In  a  more  convenient  form  as  well^by  defining 

p(^)^ 


/^T 


-  i 


(40) 


whence 


(41) 

The  self-consistency  requirement  on  the  Green  function  is  now 
contained  in  the  condition 


where 


(42) 


This  condition  should  determlnet6>(a),  and  thence  which  is 

just  i»^(o),  as  we  see  from  equation  (2ti).  It  is  more  convenient 
to  determine -Q.(s  )j  and  thence  to  flndi@/(a)  by  equation  (5^)-  In 
fact  inserting  equations  (58)  and  (59)  into  (42)  gives  a  differ¬ 
ential  equation  forJ^(a). 


(44) 


This  differential  equation  is  the  analogue  of  the  set  of  2S  coupled 

0"P 

equations^Tahir-Khell  and  ter  Haar. 

To  completely  determine  the  solution  we  require  two  boundary 
conditions.  Prom  the  definition  (56)  we  observe 

Ti(o)  =  i  (415) 

The  second  condition  is  provided  by  the  operator  identity 

TT  f)  -O  (45) 

f=-S 

Taking  the  average  of  this  equation  it  can  be  written  in  the  form 


S 

il 


TT  -A-O-O)  =  0 


(“tT) 


In  the  Appendix  it  Is  shown  that  the  solution  of  the  differential 


equation  (44)  satisfying  the  boundary  conditions  (45)  and  (47)  Is 

from  which  orl^(a)  can  be  found  by  differentiation. 


'7 


6.  The  Formal  Solution 

For  convenience  and  clarity  we  recapitulate  the  final  form  of 
the  equations,  preparatory  to  analyzing  their  low- temperature  and 
high  temperature  behavior.  Curie  temperature,  etc. 

Given  a  spin  magnitude  S,  a  temperature  T,  a  magnetic  field  H, 
and  an  exchange  Interaction  with  Fourier  components  J(K'-  the 
quasi-boson  occupation  number  <|>(<'  is 

_ I 

^  _l  (49) 

where 


This  is  an  implicit  equation  for(|)(i^  ',  involving  the  unknown  quantity 
However  ^S^is  given  in  terms  of 

$  ~  (51) 

'''  _ _  f 

(52) 

Thus  equations  (49)  -  (50)  constitute  a  set  of  coupled  equations 
which  must  be  solved  self-consistently  for 


7.  Nearest  Neighbor  Interaction 


For  simple  lattices  with  nearest  neighbor  interaction  the  for- 

9 

malism  simplifies  markedly.  In  this  case  the  exchange  J(k)  is 

(5,) 

where  d  goes  over  the  nearest  neighbors  of  a  representative  ion; 
we  assiime  all  magnetic  ions  are  crystal ographically  equivalent.  Con¬ 
sider  the  sum 


=  7  -TCk'-k  3 


K' 

—5 

which  appears  in  equation  (50)  for  E(kO*  Then 

i  fc' 


(54) 


(55) 


The  s\im  H*  is  clearly  Independent  of  S,  by  symmetry,  so 

that  it  can  be  replaced  by  I.  I.  i 

where  is  the  number  of  nearest  neighbors.  We  thus  obtain 

‘K  =  T(oiL  J'*'- (56) 


where  j(o)  = J  is  the  k  =  o  Fourier  component  of  the  exchange 
interaction.  This  equation  was  first  pointed  out  by  Mlchelene 
Bloch.  The  svim  '  Cp^tc')  is  a  function  only  of  the  temperature, 

,  and  H  (and,  of  course,  of  J  and  the  lattice  structure). 


and 


(57) 

(58) 


Thus,  for  simple  lattices  (Including  simple  cubic,  body-centered 
cubic,  and  face-centered  cubic  lattices  in  particular)  and  for 
nearest  neighbor  interaction  the  simple  spin  wave  energies  are  re¬ 
normalized  by  a  factor  depending  only  on  the  temperature  (for  zero 
field).  Independent  of  the  wave  vector  k.  Equations  (57)  and  (58), 
together  with  the  definition  of  (^(Ic)  in  terms  of  E(l<),  constitute 
a  pair  of  coupled  equations  for  the  renormalization  function  f.  For 
these  lattices  and  nearest  neighbor  interaction  the  complete  solution 
would  be  obtained  explicitly  by  solving  eqns.  (57)  and  (58)  for  f, 
thereby  obtaining  (|>(k;)  or^ln  terms  of  T,  and  H.  Then  eliminating 

^  between  this  equation  and  equation  (52)  would  give  as  a  function 

of  T  and  H. 


2  0 


B.  The  Low  Temperature  Region 

The  summation  (or  Integral)  over  R  involved  in  the  calculation 
of ^ for  H  =  o  occurs  in  simple  spin  wave  theory  and  has  been  carried 
out  by  Dyson^^^and  others  by  standard  series  expansion  procedures 
appropriate  to  low  temperatures.  Thus,  if 

E(i?)  -  2SR 


and  if  a  reduced  temperature  is  defined  by 

sA't 


4Tr<^T5  V 


(59) 


(60) 


then 


(1) 


+  -•  > 


(61) 


Similarly 


NJT 


The  constants  V  and  6j  depend  upon  the  lattice  structure  and  are 
defined  by 


V  ^  0^ 

~  33/31. 

1  Cubic. 

(63) 

^  to  r 

?or  C\j|o;C 

(64) 

V  --  f 's  , 

<0  =  4  ^or 

'  ceio'jfrf  e/  Cwlayc^ 

(65) 

To  apply  these  results  to  our  case  we  take  (compare  eqns.  (59) 


and  (53)) 


(66) 


The  remaining  equation  is  eq.  (52)  for  <^S^^in  terms  of^j 
expanding  it  in  powers  of  ^ (which  is  small  at  low  temperature), 
we  find 


(bT) 


To  disentangle  the  coupled  equations  we  first  substitute 
equation  (66)  for  R  into  equation  (62),  solving  for  f  as  a  series 
in  ‘T. 


Inserting  this  series  into  equation  (6l)  for^we  find 


(e^9) 


Finally,  inserting  this  equation  into  equation  (67)  we  find 

<S'>  5  -  ^  (I)  t’'*  -  X  ir  >  j  (£ )  r -  r' w  T'  ‘  r’''  + . . . 


.3St3/i. 


(70) 


+  (^s+i')*'  ^  (Vi) 


33+ 


Neglecting,  for  the  moment,  the  terms  \x\^ '*' 

we  see  that  the  coefficients  of  are  in  exact 

agreement  with  the  Ijyson  result.  The  term  In  is  absent,  as 

4  fl  ^ 

required  by  Dyson.  The  term  involving  is  given  by  Dyson^  •'as 


4 

where  Q,  the  ratio  the  ratio  of  this  term  to  our  t  term,  is 

-  I  [(^5  -O  ^  U 

with 


(71) 

(72) 


10,  Ur  S.  B.  .c.  ,  Fc,c. 

U  =  .31  ^  ^or  S.c^  S.cc.^  Tcc. 

(74) 

4 

Thus  our  term  corresponds  to  the  leading  term  (Q  =  l)  of  the 
Dyson  result;  that  Is,  to  the  result  which  Dyson  finds  in  the  first 
Born  approximation. 

The  terms  ^'^^,...ln  equation  (70)  seem  to 

be  spurious  results  of  the  Green  function  approximation.  For  S  =  1/2 

3  4 

they  gave  Incorrect  contributions  to  the  and  T  terms.  For  spin 
unity  the  first  spurious  contribution  appears  in  the  term,  and 

for  spin  3/2  it  appears  first  In  the  term;  thereafter  it  moves 
rapidly  to  higher  order  in  Consequently  the  spurious  terms  are 
of  consequence  In  the  low  temperature  region  only  for  spin  1/2. 

Finally,  it  is  of  Interest  to  substitute  equation  (70)  into 
equations  (b9),  (b8),  and  (bb),  to  obtain^,  f  and  the  renormali¬ 


zation  factor  R  explicitly  as  series  in  t.  We  thus  find 


^  r  *■  (75) 

-i^rrv{ii)S(i)  ^-■•■ 

-h  I--  ' 

and 

+  un  -i  ^4'V^ -  » ^S/t) 

It  will  be  recalled  that  R  is  the  "renormalization  factor";  the 

ratio  of  the  actual  quasi-particle  energies  to  the  simple  spin 

wave  energies.  At  low  temperatures  the  leading  temperature  depend- 

K/2 

ence  of  R  arises  in  the  term.  Hence  the  spin  wave  energies 


are  renormalized  by  a  factor  proportional  to  the  thermodynamic 
energy  rather  than  the  magnetization  in  this  temperature  region. 


9.  The  Curie  Temperatures,  and  the  High  Temperature  Between 

Just  below  the  Curie  temperature  (assuming  H  =  o)  the  average 
magnetization  <Cs^' is  small,  and  the  mean  niimber  of  excited  quasi- 

particles  is  large.  In  fact,  expanding  equation  (52)  in  powers  of 

-1 


i- 


we  find 


3  -  -  (78) 

Ii*urthermore ,  since  E(ic)  is  proportional  to^S^^  the  exponential  in 
the  Bose  distribution  can  be  expanded,  giving 

Multiplying  equation  (79)  by  and  replacing  by 

(from  equation  7^)  we  then  have,  in  the  limit  0 

5(sfi)  _  'hr 


i  '-irU 


(79) 


Zro)[u^<;] 


F^-O 


(80) 


where  P(-l)  denotes  the  summation 
cr/  \  -  JL  V  - 

r('i)  -  u  ^  Tlo)-T(K)  (bi) 

Ki* 

This  summation  has  been  evaluated  by  Watson;  it  has  the  values 
F^l)  ).31320(B.C,C);  (FX.c.) 


Equation  (bO)  detemlnes  the  Curie  temperature.  However  we  must 

<'S*> 

evaluate  the  limiting  value  of  the  quantity  - 


(b2) 
ust 

•f  which  appears  in 


that  equation.  Prom  equation  (57)  defining  f,  again  expanding  the 
Bose  factor. 


F(-i)-l 


-kr 


(m) 


or 


Eliminating  -V  between  this  equation  and  equation  (dO)  we  finally 
find  the  equation  determining  the  Curie  temperature. 


(«5) 


In  Table  I  we  give  the  values  of-/iT^/jas  estimated  by  Brown  and 
(3) 

Luttlnger'  -'by  extrapolation  of  the  series  expansion  for  the  sus¬ 
ceptibility  in  powers  of  1/T.  We  also  list  the  values  obtained  by 
Tahir-Klieli  and  ter  Haar^^^,  and  the  values  calculated  from  equation 
(d5).  It  will  oe  noted  that  our  values  are  higher  than  those  of 
Brown  and  Luttinger,  whereas  those  of  Tahlr-Khell  and  ter  Haar  are 
lower.  For  spin  2  the  deviation  of  our  results  from  the  Brown  and 
Luttlnger  values  is  of  the  order  of  and  the  agreement  improves 
with  increasing  spin.  In  fact  Bomb  and  Sykes  have  recently  pub- 
llshed  an  estimate  of  the  limiting  value  '^,*1-  s.  for  large  S,  for 
the  face  centered  cubic  lattice.  This  estimate  was  obtained  by  a 
painstaking  examination  of  the  systematlcs  of  the  extrapolation  of 
high  temperature  series.  They  obtain  the  value  of  6.384,  with 
which  our  value  of  6.45  agrees  within  1^. 

Domb  and  Sykes  also  give  estimates  o^Tc/j  for  S  =  •!•  and 
S  =  1  only,  for  the  face-centered  cubic  lattice;  their  values  are 
4.07  and  11.95  respectively.  These  are  lower  than  the  Brown  and 
Luttlnger  results,  and  fuirther  aggravate  the  disparity  between  our 
values  and  the  estimated  values  for  these  small  spins,  it  is 
apparent,  both  from  the  Curie  temperatuees  and  from  the  low  tempera¬ 
ture  results,  that  our  approximations  are  more  reliable  for  large 
spin,  being  particularly  bad  for  spin  Fortunately  most  cases  of 
practical  interest  are  associated  with  spins  greater  than  for  which 
the  results  appear  to  be  quite  reliable. 


»h|oj  rH  r<^lcvJ  CM  Ln|CM 
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TABLE  I 

Curie  Temperatures  for  Cubic  Lattices,  Nearest  Neighbor  Interaction 


Simple  Cubic _ |  .Body  Centered  Cubic  ),  Face  Centered  Cubic 


Brown 

Lutt, 

T-Khell 
t « Haar 

Callen 

Brown 

Lutt. 

T-Khell 
t  Haar 

Callen 

Brown 

Lutt. 

T-Khell 
t . Haar 

Callen 

1.9 

2.0 

2.7 

2.59 

2.9 

5.7 

4.2 

4.5 

5.6 

5.4 

5.5 

6.5 

7.82 

7.7 

9.1 

12.7 

11.9 

15.9 

10.6 

9.4 

11.2 

15.42 

14.4 

16.6 

24.7 

22.5 

25.5 

17.5 

15.8 

17.9 

25.17 

25.0 

26.2 

4o.O 

55.7 

4o.l 

25.8 

25.1 

26.8 

57.10 

55.5 

57.9 

58.7 

52.1 

58.5 

55.7 

51.6 

56.4 

51.19 

45.9 

51.6 

80.9 

71.4 

79.5 

For  the  face-centered  cubic  Bomb  and  Sykes ^ ^ ^glve 
Llm  I \  =  6.58  (Bomb  and  Sykes) 

VD■s(5^•o  ; 

=  5 •95  (Tahlr-Khell  and  ter  Haar) 
=  6.45  (Callen) 


^1 


It  should,  perhaps,  be  noted  that  the  Curie  temperatures  are 
quite  sensitive  to  the  decoupling  parameter  o(,  of  equation  (20). 

If  the  chosen  value  of  were  to  be  multiplied  by  (<S»>/s)*  , 

where  6  Is  any  positive  constant  no  matter  how  small,  the  Curie 
temperatures  would  become  Identical  to  those  obtained  by  Tahlr- 
Khell  and  ter  Haar,  whereas  the  low  temperatures  and  the  high 
temperature  behavior  of  the  theory  would  remain  unaltered. 

Finally,  the  high  temperature  expansion  of  the  susceptibility  % 
Is  of  Interest.  We  assme  ^S^^small,  and  maintain  only  terms  pro¬ 
portional  to  the  applied  magnetic  field.  Equation  (78)  remains 
valid,  equation  (79)  contains  the  additional  Zeeman  term,  and  the 
analogue  of  equation  (8o)  becomes 


S(sti) 

3 


nM 


[t/oj-k?)!  I  i-f- ~r 


(86) 


Similarly  the  analogue  of  equation  (83)  Is 


It  Is  convenient  to  denote 

=2: 

and  to  define  a  quantity  TC,,  by 

X  -  SZl—  ^  =  yUllLll!}  (  I  *  'X',') 


(83) 


(89) 


Then  multiplying  equation  (86)  by  [ I  4 
and  equation  (87)  by  S''J^o')X 


and  adding,  we  find 


HS4I 


X, 


~0 


3S  ■“  lf%,  (90) 

Furthermore,  expanding  the  sumKiand  in  equation  (86)  we  find 


y’  ~  'i 

In  this  summation  we  have  employed  the  identities 

F(i)  -I 


(91) 


(V 


F(z)^  1±L 
1r 


where 


T^o) 


(92) 


(95) 


The  quantity  X  can  now  be  eliminated  between  equations  (90)  and  (91 )> 
enabling  to  be  evaluated  in  a  series  in  l/T.  In  this  way  we  find 


(94) 


where  T^  is  the  Curie  temperature  of  the  molecular  field  theory: 


--  \  'irs(su) 


(95) 


The  two  leading  terms  in  this  expansion  are  in  agreement  with  the 
(^5) 

values  found' •'by  a  rigorous  expansion  of  the  susceptibility  in 
powers  of  l/T. 
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It  is  easily  corroborated  that  the  solution  of  the  differential 
equation  (44),  satisfying  the  boundary  conditions  (46)  and  (47)  is 


0^  -cjCSii^a.)  oOsCo('S,c) 

w 


(A.l) 


where 


(A. 2) 


The  evaluation  of  the  derivatives  A  CO  Ck  ^  o')  Is 

required  to  reduce  this  result  to  equation  (48).  Consider 


Is  then 


(A.5) 


where 


t" 


=  e 


d  n  *' 

Expanding  In  powers  of  y,  and  noting  that 


we  find 


(A.4) 


(A.5) 


(A.6) 


We  now  take  a  =  o  (or  y  =  l),  let  x  =  S  +  1,  and  change  variables  from 
n  to  m  =  n  -  S  and  from  p  to  r  =  S  -  p.  Then 


Similarly,  taking  a  =  o,  letting  x  =  -S  and  changing  variables  from 
ntom=n-S-l,  and  from  ptor=S+l-pwe  find 
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0) 


no/  i  $  / 


(A.8) 


Prom  equations  (A. 7)  and  (A.8)  we  note  that  the  ratio  of  the  re¬ 
levant  quantities  is 


>  ^4 1 

oO^Cj(-S,o)  (l-h$) 

oO^  6o(s+l,o') 


(A.9) 


Finally  inserting  this  ratio  into  equation  (A.l)  gives  equation  (48). 
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